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Abstract 

We present an open loop (bang-bang) scheme to control decoherence in a generic one-qubit quan- 
tum gate and implement it in a realistic simulation. The system is consistently described within the 
spin-boson model, with interactions accounting for both adiabatic and thermal decoherence. The 
external control is included from the beginning in the Hamiltonian as an independent interaction 
term. After tracing out the environment modes, reduced equations are obtained for the two-level 
system in which the effects of both decoherence and external control appear explicitly. The con- 
trols are determined exactly from the condition to eliminate decoherence, i.e. to restore unitarity. 
Numerical simulations show excellent performance and robustness of the proposed control scheme. 



PACS numbers: 03.67.Hk, 03.65.-w, 89.70.+C 



I. INTRODUCTION 



Quantum computing holds an extremely appealing practical promise and exerts a no 
less remarkable intellectual fascination. Indeed, the quantum algorithms proposed until 
now could simplify considerably the computational complexity of classically hard problems. 
While their actual implementation in robust and realistically sized quantum computers 
remains still a promise, there is general consensus in expecting further significant advances 
both in developing new quantum algorithms and the hardware they could be implemented 
in 

In the presently accepted paradigm, a quantum computation is a sequence of unitary 
transformations effected by the quantum computer upon a state in a finite dimensional 
Hilbert (i.e. a complex Euclidian) space. At the end of the computation, the system is 
left in a state which, upon measurement (read-out), would yield with high probability the 
desired result. While the read-out is non-invertible, quantum computation per se is deter- 
ministic and reversible. Thus, the realization of any quantum computer relies on maintaining 
quantum coherence in the system, for a period of time spanning at least the duration of the 
desired computation. Unfortunately, due to coupling to the environment which is not ex- 
plicitly accounted for in quantum computation, the evolution within the quantum computer 
loses its unitary character; in other words, the system decoheres 

Any quantum computation, no matter how general or convoluted, can be approximated 
to arbitrary precision using a universal set of quantum gates e.g., the discrete set of gates 
made up of the CNOT, Hadamard, phase, and it/ 8 gates |l[0]. Thus, a successful approach 
to the decoherence problem in any quantum computer requires to address and solve the 
problem at the elementary gate level. To be specific, in this paper we shall focus on the 
NOT gate, which will serve as an illustration for any one-qubit gate. The analysis of the 
two-qubit CNOT gate will be considered elsewhere. 

The NOT gate acts on a single qubit, which is the most general superposition state of a 
two-level system. Since the set containing the identity and the Pauli matrices, {J, ax, Cy, a^}, 
forms a basis in the space of 2 x 2 matrices, the coupling of a two-level system to the outside 
world can be essentially realized through the matrices az and cr^. These two couplings lead 
to different effects, known as adiabatic and thermal decoherence, respectively. 

In general, quantum decoherence is a very fast process, whose speed depends primarily 
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on the size and temperature of the computer, but may depend significantly on other factors 
as well, such as imperfect state preparation, undesired excitations, inaccuracies in manipu- 
lations of the logic gates, explicit stochasticity and noise, etc. For example, decoherence is 
induced by the spontaneous emission of a photon from a two-level atom, where the environ- 
ment is represented by the weakly-coupled electromagnetic field. Whereas the evolution of 
the two-level system and its environment considered together as a large (closed) quantum 
system is in principle unitary, the subdynamics of the two-level system generally is not. 
Clearly, this behavior represents a major stumbling block on the path to realizing a usable 
quantum computer at room temperature. 

Several approaches have been proposed to eliminate or mitigate the undesirable effects of 
decoherence in open quantum systems, including open loop (quantum bang-bang) control 
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and quantum feedback 



ree subspaces (DFS) ll|, , quantum error correction jisl . lJ | 
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The open loop (quantum bang-bang) control has been pioneered by the seminal work of 
Lloyd and Viola (see Refs. 0,y,[lO| and references therein). Vitali and Tombesi [3] have also 
considered applying a sequence of frequent parity kicks, as well as an appropriate stochastic 
modulation, in order to control decoherence. Essentially this type of approach relies on 
applying rapid control pulses to the system in order to dynamically cancel the system- 
environment interaction. It has been shown that decoherence can be effectively suppressed 
if the pulse rate is much higher than the frequency corresponding to the correlation time of 
the environment. Thus, this open- loop control technique "freezes" decoherence in a manner 
analogous to the quantum Zeno effect. 

DFS are special subspaces of the Hilbert space, which are inherently immune to decoher- 
ence, due to specific symmetries in the system-environment interaction. This concept was 



first introduced by Zanardi and Rasetti 
subsystems [14J. It has also been shown 
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, and has been recently generalized to noiseless 
Q that bang-bang control can be used to cre- 



ize a quantum state has been reported for photon singlets .IE], for 



ate the conditions required for large DFS. The successful use of decoherence free subspaces 
(DFS) to passively stabi 
a pair of trapped ions [20| , and more recently for a noiseless subsystem via liquid state NMR 
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Quantum error correction is historically the first method proposed to counter decoher- 
ence, and also the only one to date for which an accuracy threshold theorem has been 
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proved [2^. It is essentially an active form of quantum feedback control implemented on a 
system encoding redundant information. Valid states of the system are restricted to special 
Hilbert subspaces (codes) so that any loss of information due to decoherence (errors) can 
be diagnosed and reversed. The system is closely monitored to observe errors, which can 
be corrected later by conditionally performing appropriate feedback operations. It has been 
shown that fault-tolerant recovery which allows for errors in the correcting feedback, if the 
error rate is below a critical value. 

The efficient implementation of both DFS and quantum error correction techniques re- 
quires a large amount of ancillary quantum memory. Thus quantum feedback (closed-loop) 
and open-loop control strategies have the advantage that in general they do not require any 
ancillary quantum resources. 

Quantum feedback can also be used to stabilize a quantum state against decoherence, as 
first shown by Tombesi et al. in a series of papers [l5|. For example, in a scheme recently 
proposed by Wang et al. |3l for a two-level atom, photocurrent feedback from homodyne 
detection of spontaneous emission is used to alter the atom-environment dynamics. The 
quantum state of the atom is driven to and then remains in a pure target state, even though 
the atom is still spontaneously emitting. Of course, the feedback delay time has to be shorter 
than the time scale corresponding to the spontaneous emission rate, in order to successfully 
maintain the atom in the target state. We note that this quantum feedback scheme cannot 
be used to maintain an unknown state such as that obtained from a quantum computation, 
it can only drive the system to a known target state. 

In this paper, an open-loop control scheme is proposed and simulated, to eliminate the 
effects of decoherence in a two-level system. While the general philosophy is following in 
the steps of Lloyd and Viola's seminal results 0,0, 3, our implementation has several new 
aspects. More precisely, we present the first implementation of a quantum control scheme in 
which: (i) decoherence and control are taken to act simultaneously within a realistic model, 
which allows one to deal with thermal decoherence; (ii) the required control is directly related 
to and calculated from the decoherence effects, which presents the practical advantage of 
maintaining the frequency and amplitude of the required controls at manageable levels; and 
(iii) the effect of imperfect control pulses is assessed. We are aware that determining the 
decoherence behavior of a quantum system can be difficult in practice, but if we know it 
then our control scheme can take advantage of this, to tailor more efficient controls. 
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The remaining of the paper is organized as follows: In section 2 we present the spin-boson 
model for a two-level system in contact with the environment and an external control. In 
Section 3 we calculate the reduced density matrices for the two level system after tracing 
out the environment modes, in the case of adiabatic and thermal decoherence. The control 
strategy is implemented in Section 4 and the results are discussed in Section 5. 



II. SPIN-BOSON MODEL 



The exact quantum mechanical description of the system starts from the complete Hamil- 
tonian that accounts for the two-dimensional system and the environment which, in principle, 
is the rest of the Universe. In different situations the environment predominantly manifests 
itself in different ways, e.g., electromagnetic modes, acoustic modes, etc. Thus, in many 
instances, a reasonable and generic model is represented as a system of noninteracting bo- 
son modes. The two-level system interacting with these modes is known as the spin-boson 

n 

model p|. We write the complete Hamiltonian in the form: 

H = Hs + He + Hi + H, (1) 

where 

2 

Hs = J2 PiiEi (2) 

i=l 

and 

oo 

H^ = h^ ujQqa\aq (3) 

q=l 

represent the two-level system and the boson bath. Here Pa denote the projection operators 
Pa = and Ei are the corresponding energies. The terms 

oo 
<?=! 

and 

He = — hQpV(t) a xcos{ujot) (5) 

describe the interaction of the two level system with the environment and the external 
control, V{t), respectively. The interaction with the environment is parametrized by e 
(whose magnitude indicates the strength of the coupling) and can result in a phase flip 
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or a bit flip (a bit-phase flip can also be obtained either by combining the two effects or 
by including a ay term). However, from a formal viewpoint, the effects of ax and ay are 
essentially similar and we shall not consider the latter. On the other hand, the control 
which is applied as a fast laser pulse, is assumed to act strongly only upon the spin, without 
signiflcantly affecting the environment, either directly or indirectly. 



While the spin-boson model in various guises has been thoroughly studied 
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26l |. the consistent inclusion of the control term. He-, in the Hamiltonian, and the general 



derivation of the corresponding reduced description are new. 

The dynamics of the two-level system considered in this paper is expressed in terms of 
parameters normalized to the Rabi frequency, Vtp, namely the time evolution is described 
in terms of dimensionless Rabi time units, t = VLpt, and all the frequencies, namely ujq = 
ujoq/^Fi ^ = ^o/^Fi Qq = ^q/^F and Co'12 = Qu/^F are also renormalized with respect to 
the Rabi frequency. The frequency Q12 is the transition frequency between the energy levels 
of the two-level system. 

The state of the system at r = is described by the density matrix: 



p(s,e,0) =p(s,0)®p(e,0) (6) 

where 

pie,0) = l[pq{e,0) (7) 

p,(e,0) = [1 -e-^]-i5^e-^"«|n,)K| (8) 

riq 

and 

2 

p{s,0) = J2p^MP^r (9) 

In general, the decoherence experienced by a two-level system has two distinct timescales, 
corresponding to the adiabatic and thermal regimes separately. 

Adiabatic decoherence is responsible for the decay of transverse (x, y) polarization (the 
off-diagonal density matrix elements). It acts on a relatively short timescale, such that 
a superposition state will typically decay to a mixture of states. For larger, macroscopic 
sub-systems, adiabatic decoherence ensures that quantum superpositions of distinct states 
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will not be observed. The mixture of states generated by adiabatic decoherence will then 
continue to relax to a stationary value over a longer thermal timescale. 

Thermal decoherence generally causes a change of all density matrix elements, and leads 
to the exponential decrease of the excited state population of the two-level system, due 
to energy exchanges between the two- level system and its environment. It is commonly 
neglected in discussing decoherence of quantum systems since its timescale is longer than 
the timescale for adiabatic decoherence, which is deemed responsible for eliminating most 
quantum effects [27|. We note that often thermal decoherence is refered to as (quantum) 
dissipation, and adiabatic decoherence is simply refered to as decoherence. 

To simplify the calculations and render them more transparent, it is convenient to write 
the evolution of the system in the interaction representation, and consider separately the 
adiabatic {a = 0, /3 = 1) and thermal [a = 1, (3 = 0) decoherence regimes. A combined 
treatment of these regimes is straightforward, though computationally cumbersome. The 
unitary evolution operator in the interaction representation, Uj, is related to the unitary 
evolution operator in the Schrodinger representation, U, by 

Ui{s, e, r) = e^^o"/^f/(s, e, r) (10) 

where Hq = Hg + Hg. 

Thus the evolution of the density matrix of the system and environment is 

where the interaction Hamiltonian contains both the interaction with the environment and 
with the applied control: 

Hj = Hj, + Hn, (12) 



Hic = ~V{T)a,, (13) 



-eha,j:^{g;age-^^''^ + g^ale'^^^n a = 0, P = 1 

Hr = \ (14) 

-eh[^^ E,^?,Xe-*(---')- + E,^7</4e^^-"""'^")] a = 1, /3 = 

These interaction-picture Hamiltonians have been calculated by using the rotating wave 
approximation, and assuming a zero detuning 5 = uj\2 — = 0. 



The evolution equation for p admits the formal solution 

p(s, e, r) = Ui{s, e, r)p(s, e, 0)[/|(s, e, r) (15) 

where the evolution operator satisfies the equation 

-j^ = ^^i^^c + Hji\Ui{s, e, r) (16) 

given the initial condition Uj{s, e, 0) = 1. The formal solution of this equation is 28^ . 

Uj{s, e, t) = r[exp{-^ dr'iHUr') + Hnir'))}], (17) 

where T[ ] represents a time-ordering operator. 

Since the various terms in the interaction Hamiltonian Hj do not commute, the complete 
evolution operator, Uj{s,e^), cannot be calculated exactly, except in special cases. For 
instance, Viola and Lloyd jSj have shown that if the interaction with the environment Hn 
contains only cr^, the evolution operator Uj{s,e,T) can be substantially simplified by using 
vr-pulses for the control. However this type of control would not allow the same simplification 
if the interaction with the environment Hn contains ax or ay. 



The general Baker-Hausdorff theorem [291] can be used to expand the evolution operator 
into an infinite product of exponentials, 

Ui{s, e, r) = e^^^o ^tHi^it) ^ dtHnit) ^ j; dtj^ M[HUt),Huit')\ ^ _ _ _ (^^g) 

We assume that the effect of the control and interaction Hamiltonians is relatively small 
at all times. Neglecting the commutators, the evolution operator can be written as a first 
order approximation in the magnitude of the control pulses, V{t), and the coupling strength 
parameter, e, of the system-environment interaction. 



f//(s,e,r) 



exp{-{ £ dtHjcit)} X exp{-| dtHnit)} « = 0, /5 = 1 

exp{-i £ dt{Hi,{t) + HnSt))} x expj-f £ dtHn^{t)} a = 1,(3 = 



(19) 

where Hn^ and Hn^ correspond to the terms in the system-environment interaction Hamil- 
tonian Hji, proportional to a^ and ay respectively. 

Performing the time integrals in the exponents, we obtain 

exp{?cr^./(r)} X exp{— 2cr2Q_(r)} a = 0, (3 = 1 

Ujis,e,r)={ (20) 

exp{io-^/(r) - axQ-ir)} x exp{iayQ+{T)} a = 1, (3 = 
8 



where 



/(r) = \ / rfrV(r') (21) 







Q-{r)='-Y,(^M,a\-Mla,) (22) 
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Q+(r) = i^(M,at + M;a,) (23) 



2 

1 



and 



M,(r) (24) 

_ Qi{u^l2-^^)r\ ct = l, /? = 

III. EVOLUTION OF THE REDUCED DENSITY MATRIX 

To evaluate the various effects of decoherence we have to calculate the reduced density 
matrix of the two-level system, by tracing out the environment modes: 

2 

Pi,{s,T)=Tr,{Y,^\UiPkip{e,m\\j)Pkm} (25) 

k,l=l 

for i,i = 1,2. 

A. Adiabatic case (a = 0, /? = 1) 

Upon expansion of the exponentials in Eq. (j2(Jj) . we obtain the following approximate 
expressions for the evolution operator and its adjoint: 

Ui{s, e, t) = Eo + iE^a^, + Eyay + E^a^ (26) 

and 

Uj (s, e, t) = Eq- iE^a^ - Eyay - E^a^, (27) 

where 

Eq = cosh(2g_(r))cos(/(r)) (28) 

E^ = cosh(2g_(r))sin(/(r)) (29) 

Ey = -sinh(2Q_(r))sin(/(r)) (30) 

E, = -sinh(2g_(r))cos(/(r)), (31) 
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with /(r) and Q_(r) defined by Eqs. (PT|)-(P^. 

Using these expressions in Eq. fj^ . we get the reduced density matrix of the two-level 
system: 



Pii = Pii(0)cos2/ + p22(0)sin2/-z[pi2(0) -p2i(0)]e-3"''cos/sinJ, (32) 

P22 = P22(0)cos2/ + pn(0)sin2/ + z[pi2(0) -p2i(0)]e-f«^cosJsin/, (33) 

Pi2 = Pi2(0)e-^'»^cos2/ + p2i(0)e-^''"'sin2/ + i(p22(0) -pii(0))cosJsin/, (34) 

P2i = P2i(0)e-^"^cos2/ + pi2(0)e-3'"'sin2/-^(p22(0) -pii(0))cosJsin/, (35) 

where / = /(t), and 

Qad ■= gad{r) = 7 / duG{u){l - cosut) coth ^ (36) 
Jo ^ 

is the decoherence function obtained by Palma 2^. The dimensionless constant 7 depends 

on the dipole moment of the two-level system and on the Rabi frequency. The strength 

of the decoherence rate experienced by the two- level system is determined by 7, ie., the 

decoherence rate is weak for 7 << 1. 

The function GiiS) is the spectral function 



where n is the dimensionality of the system (usually n = 3), and Uc is the (usually very 
large) cutoff frequency, which ensures the convergence of the improper integral. In the 
case of phonon modes this cutoff frequency is the Debye frequency. We note that the only 
dependence on the two-level system is contained in 7. 

The expressions for the reduced density matrix elements p^ = 1,2) retain 

the correct limit behaviors, and the trace of the density matrix is always equal to one. 
Indeed, in the limit of zero control we recover Palma's result [2^, whereby only the non- 
diagonal elements are affected by the decoherence process, while the diagonal elements 
remain unchanged. When both the decoherence and the control are equal to zero, we 
recover the ideal situation in which the state remains unchanged. 
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B. Thermal decoherence (a = = 0) 



Proceeding as before, we use the approximation in Eq. for a = 1 and /3 = 0, to 

obtain the following expressions for the evolution operator and its adjoint in the thermal 
decoherence case: 

Ui{s, e, t) = Eo + iE^a^ + EyCy + E^a^ (38) 

and 

Uj{s, e, t) = Eq- iE^a^ - iEyay + iE^^a^ (39) 

where 

Eq = cos(/(r)+«g_(r))cos(g+(r)) (40) 

E, = sin(/(r) + tQ4r)) cos(g+(r)) (41) 

Ey = cos(/(r)+zQ_(r))sin(g+(r)) (42) 

E, = sin(/(r) + zQ_(r)) sin(Q+(r)), (43) 

with J(r), g_(r) and Q+{t) defined by Eqs.(j2I])-(|2Sl)- To calculate Uj{s, e, r), we have used 
the fact that Q~{t) and <5+(t) commute in the first order approximation for the coupling 
strength parameter, e. 

The environment modes are then traced out to obtain the elements of the reduced density 
matrix: 

Pii = ^[1 + (Pii(O) - P22(0))e-2^"' cos(2/) - z(pi2(0) - P2i(0))e-^*'' sin(2/)] (44) 
Pi2 = [i?e{pi2(0)} + /m{pi2(0)}cos(2/)]e-^--^[(pn(0)-p22(0))e-2^-sin(2J)] (45) 
P21 = [i?e{p2i(0)} + /m{p2i(0)}cos(2/)]e-^'- + ^[(pn(0)-p22(0))e-2^-sin(2/)] (46) 
P22 = ^[1 - (Pii(O) - P22(0))e-2^"^ ^^g(2/) + z(pi2(0) - P2i(0))e-^*'^ sin(2/)] (47) 

where the decoherence function, gth, is given by 

9tH := 9th{T) = 7 ^^ l-^^^[(^i2- ^)^] ^3 eoth(/?oc^/2) exp(-.;M) (48) 
In the limit of zero control, the matrix elements read: 
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Pl2 



^(l + e-2^-(pn(0)-p22(0))) 
^12(0)6-^^*" 



(49) 
(50) 
(51) 
(52) 



P21 



P2i(0)e- 



P22 




(Pll(0)-P22(0))) 



IV. CONTROL STRATEGY 

Our control strategy is based on the following idea: by equating the elements of the 
reduced density matrix in the presence of decoherence (adiabatic or thermal) and (unknown) 
control (Eqs. (jH^ - ljH^j) or (03)" (03) 5 respectively) with the elements of the density matrix 
undergoing a unitary evolution, we can, in principle, determine the control that eliminates 
the effect of the decoherence and momentarily restores unitarity. In our case, this simply 
means that the matrix elements should be restored to their initial values, i.e. the unitary 
evolution is the identity operator. Indeed, since the primary goal here is to show robust 
elimination of decoherence, we consider that the NOT transformation has already been 
effected and the state is now waiting to be involved in the next operation, as specified by 
the ongoing quantum algorithm. Since this "waiting time" may be much longer that the time 
needed to realize the "vr-pulse" of the NOT gate, we simply ignore the latter and concentrate 
on the former. Inclusion of the "vr-pulse" due to a NOT gate has been actually implemented 
and does not lead to any significant change other than an unnecessay complication of the 
formalism. This separation of the quantum gate operations from the decoherence control 
operations allows us to treat any one-qubit gate. 

To numerically implement the control strategy, the real and imaginary parts of the four 
complex elements of the density matrix calculated in Section IIIII are used in a prescribed 
order (see below). While this particular order is not essential, we shall define and use it 
consistently here, to make it easier to follow our control strategy. Indeed, since a single 
control cannot realize the desired effect (namely the instantaneous restoration of the ideal 
behavior) for all eight elements at once, we have to adjust these elements in turn. The 
control cycle is made up of 8 control steps, since there are 8 equations to solve for the (real 
and imaginary) components of the 2x2 complex density matrix. This leads to a sequence 
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of eight real transcendental equations to be solved in turn. 

Ideally, the density matrix of a 2 x 2 system is usually completely described by three 
independent quantities. It may happen that, even in the non-ideal case, for particular 
Hamiltonians, certain symmetries impose certain relations between the matrix elements. 
However, to maintain a general and systematic character of the approach, we prefer to treat 
them all as independent variables, at the expense of a slight lengthening of the control 
process. Since this process eventually stabilizes into a rather short cycle, we do not consider 
this to be a serious drawback. 

The algorithm is apphed identically for either the adiabatic or the thermal case; thus, 
we describe it for a generic decoherence function, denoted g. In the first cycle of eight time 
steps, we start by considering pii^(l) := pii^(5f(l), J(l)) after one time step. To determine 
the control pulse required to set piifl(l) equal to its unitary value after applying a NOT 
gate, we find a control value 7(1) which solves the equation = Piiii(O). 

Finding the appropriate solutions of the above transcendental equation requires a certain 
care. Indeed, if there are no solutions in a given interval, it is impossible to restore the unitary 
behavior. On the other hand, if there are multiple solutions, it becomes very difficult to 
guarantee that the same control is used for every cycle. To obtain a unique and consistent 
solution, we had to choose carefully a specific interval. Of course, this caveat applies for the 
solutions of the other equations as well (see below) . 

Denoting the solution of this equation by /(I), we reset piiij(l) exactly to its original 
value by using a control pulse /(I). After the first time step, all the other elements will have 
suffered the effect of decoherence, and the effect of the control pulse, /(I): 



(53) 
(54) 
(55) 
(56) 
(57) 
(58) 
(59) 
(60) 




Pi2r{1) = Pi2ii(5(l),/(1)) 
Pi2/(1) = pm(^(l),/(l)) 




P21/(l) = P217(^(l),/(1)) 
P22r(1) = p22R{g{l)J{l)) 
P22/(l) = P227(5(1),/(1)) 
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By definition, the restoration of unitarity executed by the control / is exact. 

The control strategy advances through the next seven time steps in a similar fashion. 
The required control pulses are determined by solving each of the other seven density matrix 
equations in the order presented above. 

For example, at the second time step, we consider the equation for pii/(2) = 
Pii/(5'(2), /(I) + -^(2)). To determine the control pulse required to set piu{2) equal to its 
desired value, we find the control 7(2) which solves the equation piii{2) — piii{0). After 
applying this control, the values of the matrix elements after the second time step are given 
by: 



Piiii(2) 


= puR{9{i),m) 


(61) 


Pii/(2) 


= pn/(^?(2),/(l)+/(2))=pn/(0) 


(62) 


Pl2il(2) 


= Pi2k(^(2),/(1) + /(2)) 


(63) 


Pl2/(2) 


= pi27(^(2),/(l) + 7(2)) 


(64) 


P2ir(2) 


= P2i/j('?(2),/(l) + /(2)) 


(65) 


P21/(2) 


= P2i7(5(2),/(l) + /(2)) 


(66) 


P22il(2) 


= P22ij(^(2),/(l) + /(2)) 


(67) 


P22i-(2) 


= p22/(5(2),/(l) + /(2)) 


(68) 



Note that the decoherence function used to calculate Piir(2) is different from those for 
all the other matrix elements, since pii/j(l) was reset to piii?(0) after the first time step. In 
general, resetting any one of the matrix elements to its corresponding original value implies 
that the cumulative effect of the control pulses applied since the last correction (eight time 
steps in the past) cancels out the effect of decoherence since the last correction. 

To complete the first cycle we proceed to solve the equations for the next six matrix 
elements, to determine the required control pulses. By the end of the first cycle of eight 
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time steps each of the matrix elements depends on different decoherence and control values: 



Piii?(8) 


= piiRi9i7), /(2) + /(3) + /(4) + /(5) + /(6) + 1(7) + 1(8)) 


(69) 


Pii/(8) 


= Pii7(^(6), 7(3) + 7(4) + 7(5) + 7(6) + 7(7) + 7(8)) 


(70) 


Pi2_r(8) 


= Pi2r(^(5), 7(4) + 7(5) + 7(6) + 7(7) + 7(8)) 


(71) 


Pl2l[o) 


- Pi2i[g[4:), + + + I [b)) 


[72) 


P2ir(8) 


= P2ir(^(3),7(6) + 7(7)+7(8)) 


(73) 


P21/(8) 


= P2i/(^7(2),7(7) + 7(8)) 


(74) 


P22il(8) 


= P22il(^(l),/(8)) 


(75) 


P22/(8) 


= P227(^(8), 7(1) + 7(2) + 7(3) + 7(4) + 7(5) + 7(6) + 7(7) + 7(8)) 


= P22/(0I76) 



a pattern which will recur in every cycle of eight time steps. 

After one time step in the second cycle of eight time steps, we are back at the starting 
point. To determine the control pulse required to set piii?(9) equal to its initial value, we 
find a control value 7(9) which solves the equation piiji{g(8), I{2) + 7(3) + 7(4) + 7(5) + 
7(6) + 7(7) + 7(8) + 7(9))) = punigil), /(I)) = Piiij(O). 

In the following time step, set pii/(10) = pui{g{8), I{3) + 7(4) + 7(5) + 7(6) + 7(7) + 
7(8) + 7(9) + 7(10)), and so on. 

This procedure is repeated as long as the quantum state of the two-level system has to be 
maintained, i.e. the "waiting time" mentioned previously until the next logic gate is applied. 
We note two important things for the apphcability of the scheme. First, the knowledge of 
the decoherence function is needed only for a finite period of time (in the example above 
only eight time steps). Second, after initial transients, controls will stabilize and the whole 
control cycle will repeat itself periodically. In other words, for each sequence the controls 
can be calculated off-line and applied in the required order. This behavior has indeed been 
observed: the values of the control pulses were rapidly stabihzed after the first cycle of 8 
time steps, as shown in Figures 1 and 2. The stabilization time depends on initial conditions, 
and the magnitude of the decoherence. 
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V. DISCUSSION OF THE RESULTS 



The graphs shown in Figs. 1 and 2 illustrate the typical evolution of the matrix elements 
while applying this control strategy in both adiabatic and thermal decoherence situations. 

As expected, the relative size of the time step between control pulses determines the 
amount of deviation of the matrix elements from their unitary values. Also, the size and fre- 
quency of the control needed to restore perfectly the ideal situation depends on the strength 
of the decoherence. For 7 = 1 (extremely strong decoherence) , the frequency of the control 
can be decreased only to about twice the Rabi frequency. If we decrease the frequency even 
more, we cannot restore exact unitarity, at least not by using this scheme. Since the control 
algorithm described here is based on periodic perfect restoration of the ideal situation, we 
shall not discuss this case any further. However, imperfect restoration of unitarity is likely 
to happen and is not at all hopeless: this situation will be analyzed in future work. As 
the strength of the decoherence, 7, decreases (strong to medium decoherence), we can also 
decrease the frequency of the control pulses while keeping the same amount of deviation of 
the matrix elements from their unitary values. Of course, higher frequencies result in better 
restorations of unitarity. 

The same observations apply for thermal decoherence. Since the latter affects all elements 
of the density matrix, it is reasonable to expect that the same amount of control will be less 
efficient here than in the adiabatic situation. This is indeed the illustrated in Fig. 

2. 

If the controls are perturbed by noise, the control is still very efective. It is important 
to note that the amplitude of the noise is calibrated with respect to the state and not with 
respect to the amphtude of the control needed in the noiseless situation, which in some cases 
is extremely small. 

In conclusion, we presented an open loop control scheme for quantum gates which is: (i) 
realistic - the control and decoherence act simultaneously, and the model is able to deal with 
adiabatic as well as thermal decoherence; (ii) relatively simple to calculate, (iii) efficient to 
implement, and (iv) robust to noise both in the decoherence functions as well as in the 
control pulses. 

The results presented in this paper depend, among other factors, on the model Hamilto- 
nian, and the approximation used in the tracing out procedure. However, modifications of 
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FIG. 1: The unitary, adiabatically decohered, and controlled evolution of the Im{pi2} element of 
the density matrix for the initial state > "t""^!^ '^^^ dimensionless decoherence rate is set 
to strong, 7 = 1, and the other two parameters characterizing these plots are: the time between 
control pulses, T (scaled in terms of the Rabi frequency), and the standard deviation. A/, of the 
control pulses after adding normally-distributed noise, (a) T = 0.5, AI = 0; (b) T = 0.5, AI = 0.1; 
(c) r = 2, A/ = 0; (d) T = 2, A/ = 0.1. 

either would not alter the main features of the proposed control strategy. 

While our control scheme shares several common features with the quantum feedback 
scheme proposed by Wang et al. [l6[ and the quantum bang-bang control, there are also 
significant differences between our scheme and these approaches, which we briefly discuss 
here. 

The quantum feedback eliminates specifically the decoherence effects of spontaneous emis- 
sion, whereas the present scheme is more general. Both the feedback scheme and the ap- 
proach proposed here aim to restore (at least approximately) the unitarity i.e., maintain the 
information contained in the quantum state. However our scheme assumes prior knowledge 
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FIG. 2: The unitary, thermally decohered, and controlled evolution of the Im{pi2} element of the 
density matrix for the initial state > +"^1^ ^- '^^^ dimensionless decoherence rate is set 

to strong, 7 = 1, and the other two parameters characterizing these plots are: the time between 
control pulses, T (scaled in terms of the Rabi frequency), the decoherence rate, 7 (dimensionless), 
and the standard deviation. A/, of the control pulses after adding normally-distributed noise, (a) 
T=0.25, AI = 0; (b) T = 0.25, AI = 0.1; (c) T = 1, A/ = 0; (d) T = 1, A/ = 0.1. 

of the initial state in order to maintain it unchanged, whereas the quantum feedback scheme 
drives the state to a target state regardless of initial conditions. Finally, our scheme does 
not require any measurements on the two- level system or the environment, and does not 
involve any feedback. 

Si.na...y .0 ,ua,.u. bang-bang control .cben,es BflaQ. the control proposed bere 

lies on a series of relatively fast control pulses to dynamically cancel the system-environment 
interaction. However, in our approach, we tailor the control according to the specific na- 
ture of decoherence present in a physical system. This has certain advantages, namely: (i) 
controls are applied only at the required frequency and magnitude, without overtaxing the 
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available resources; (ii) the limits of applicability are also clear, beyond which a specific 
open loop scheme would not be able to restore exact unitarity; (iii) since the environment is 
not supposed to change significantly during the quantum computation, the control sequence 
stabilizes very quickly to a set of periodic controls, as shown in Figs. 1 and 2 (this implies 
that in most situations of interest the controls can be computed off-line and then simply and 
efficiently applied in the prescribed sequence); and (iv) robustness to noise appears quite 
strong, both in the adiabatic and thermal cases. 

The main difference between our scheme and the bang-bang control scheme proposed by 
Viola and Lloyd 8] is that, on the one hand we can handle the more difficult case involving 
thermal decoherence, and on the other hand the control pulses are tailored to the task 
instead of being generic. Indeed, a cycle of two predetermined 7r/2-pulses is used in Viola 
and Lloyd's scheme for adiabatic decoherence. The evolution of these pulses is disjoint from 
the evolution of the system and the environment, in order to ensure that the cycle evolution 
can be simplified to an expression containing only o"^. However, the control proposed in 
Ref. yi cannot be used as such in the case of thermal decoherence, which has a system-bath 
interaction Hamiltonian proportional to cr^^, even if the control Hamiltonian is taken to be a 
linear combination of a^^ and Oy. On the contrary, in our scheme the evolution of the pulses 
is not disjoint from the evolution of the system and the environment, which allows us to 
handle a general combination of thermal and adiabatic decoherence by employing the BCH 
approximation. Moreover our controls do not have a 'priori assigned values. In general, the 
stabihzed pulses have a much smaller amplitude than the 7r/2-pulses, for the same pulse 
duration, and they are applied at a slower rate. In principle, this improves the efficiency of 
the control. 

The limitation of our scheme is its reliance on knowledge of the target state that we 
want to maintain, e.g., the initial values of the reduced density matrix of the two-level 
system. Since the target state changes during the course of the quantum computation, the 
applicability of this control scheme is presently limited to short calculations. We intend 
to explore this aspect together with different Hamiltonians, approximations, and control 
schemes in the future. 
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